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ABSTRACT 

Protostellar systems, ranging from low-luminosity classical T Tauri and Herbig Ac stars to high- 
luminosity Hcrbig Be stars, exhibit a near-infrared (NIR) excess in their spectra that is dominated 
by a bump in the monochromatic luminosity with a peak near 3 ^m. The bump can be approximated 
by a thermal emission component of temperature '~ 1500 K that is of the order of the sublimation 
temperature of interstellar dust grains. In the currently popular "puffed-up rim" scenario, the bump 
represents stellar radiation that propagates through the optically thin inner region of the surrounding 
accretion disk and is absorbed and reemitted by the dust that resides just beyond the dust subli- 
mation radius rg^b- However, this model cannot account for the strongest bumps measured in these 
sources, and it predicts a pronounced secondary bounce in the interferometric visibility curve that is 
not observed. In this paper we present an alternative interpretation, which attributes the bump to 
reemission of stellar radiation by dust that is uplifted from the disk by a centrifugally driven wind. 
Winds of this type are a leading candidate for the origin of the strong outflows associated with proto- 
stars, and there is observational evidence for disk winds originating on scales Tgub- Using a newly 
constructed Monte Carlo radiative transfer code and focusing on low-luminosity sources, we show that 
this model can account for the NIR excess emission even in bright Herbig Ae stars such as AB Auriga 
and MWC 275, and that it successfully reproduces the basic features of the visibilities measured in 
these protostars. We argue that a robust dusty outflow in these sources could be self-limiting — 
through shielding of the stellar FUV photons — to a relatively narrow launching region between rgub 
and ~ 2rsub- We also suggest that the NIR and scattered-light variability exhibited by a source like 
MWC 275 can be attributed in this picture to the uplifting of dust clouds from the disk. 
Subject headings: circumstellar matter — ISM: jets and outflows — magnetohydrodynamics (MHD) 
— protoplanetary disks — radiative transfer — stars: individual (AB Auriga, 
MWC 275) 



1. INTRODUCTION 

Protostars ranging from low to high masses are un- 
derstood to form in molecular cloud cores that un- 
dergo gravitational collapse. In general, the gas in the 
cloud possesses angular momentum, so a core under- 
going collapse eventually encounters a centrifugal bar- 
rier and settles into a rotationally supported disk. It is 
believed that most of the mass that is assembled into 
young stars reaches the stellar surface through the as- 
sociated disk (e.g., iCalvet et al.l l20dol ). The existence 
and properties of protostellar disks have been inferred 
from a var iety of spectroscopic and imaging observa- 
tions (e.g., IMillan-Gabet etahl I2007|) . whereas the ar- 
rival of disk material onto the stellar surface (where it 
liberates its gravitational potential energy) has been de- 
duced from measurements of excess (over the intrinsic 
stellar radiatic)n) op t ical/UV continuum em ission (e.g., 
iHartigan et"allll995t [Corcoran fc Ravlll998[ ). Protostel- 
lar spectra also exhibit a near-infrared (NIR) excess, 
which can in part be attributed to the emission from 
the disk surface (representing the intrinsic disk radiation 
as well as radiation originating in the star that is ab- 
sorbed and reprocessed by the disk). However, in many 
cases involving both intermediate-mass Her big Ae/Be 
(HAeBe) stars (e.g., iHillenbrand et aH 1992f) and low- 
mass classical T Tauri (cTT) stars (e.g.. IMuzerolle et al.l 
[200l . the NIR excess is dominated by a prominent 



"bump" in the monochromatic luminosity that peaks 
near 3 fim and can be approximated by a thermal emis- 
sio n component of tempe r ature ^ 1500 K. As reviewed 
in IDullemond fc Monmerl (j2010l hereafter DM10), this 
bump cannot be reproduced by disk models that are ev- 
erywhere optically thick to the stellar radiation. 

Protostellar systems generally also exhibit strong out- 
flows that are evidently powered by the accretion pro- 
cess. It is commonly believed that these outflows 
are driven centri fugally by a large-scale, ordered mag- 
netic field (e.g., iPudritz fc NormanI |1983() . although it 
has not yet been conclusively establis hed whether this 
field originates in the protostar (e.g ., iShu et all 120001 : 
iRomanova et al.ll2009HCranmedl2009l ). in which case the 
outflow must emana te close to the protostel l ar surface, 
or in the disk (e.g., iKonigl fc Pildnti [2000l: IRav et al.l 
l2007t) . in which case the launching region can lie further 
out. If the wind-driving region extends beyond the dust 
sublimation radius Tgub (the disk radius where the equi- 
librium temperature of dust grains that are exposed to 
the protostellar radiation field is equal to the sublimation 
temperature of the grains) then disk outflows of this type 
offer a promising means of accounting for the observed 
~ 3 /im bump. This is because centrifugally driven winds 
can efficiently uplift dust from the disk (through gas- 
grain collisions), which gives rise to a vertically strati- 
fied dusty gas distributio n ab ove the disk at radii > rsub 
(|Safiedll993at see Section [Z3|) . The dusty outflow can in 
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principle intercept a significant fraction of the continuum 
radiation emitted by the protostar and the inner disk, 
which is, in turn, reprocessed to infrared wavelengths. 
The resulting spectrum exhibits a bump near 3 fim on 
account of the fact that the dust sublimation temper- 
ature is Tsub ~ 1500 K, which implies a cutoff in the 
dust emission below ^ 2 ^m, and that grains have a rela- 
tively high emiss ivity in the few-micron range. This was 
demonstrated bv lKonigl fc Karti^ (1 19941) using th e semi - 
analytic disk-wind model of iBlandford &: Pavnd ()1982l ) 
and a radiative transfer code based on the "lambda op- 
erator " perturbation technique . 

The iKonigl fc Kartid ()1994D work was carried out in 
the context of active galactic nuclei (AGNs), where a sim- 
ilar bum p had been iden tified as a common NIR spectral 
feature ( jEdelson fc Mal kan 1986), and where the con- 
nection between the wavelength of the bump and the 
expect ed sublimation t emperature of dust had been first 
noted (H arvaini i [19901) . However, in view of the basic 
physical principles that underlie this picture, it is clear 
that the "dusty hydromagnetic disk wind" mechanism for 
generating the NIR spectral bump is quite general and 
does not depend on the specific astrophysical context in 
which it is realized. Indeed, at about the same time that 
the AGN application was worked out, P. Safier used a 
simplified disk-wind model to demonstrate that this sce- 
nario could also account for the NIR spectra of proto- 
stars. However, his results were only presente d in a few 
conference proceedings (see, e.g.. lKon igr'1996') and have 
not circulated through the general protostcllar research 
communityQ A few years later, another model was pro- 
posed to explain the NIR spectra in protostars, nam ely 
the "puffed-up inner rirn " scenario (jNatta et al.ll2001t see 
also iTuthill et al.l l2001| ) , and this has become the most 
widely discussed explanation of the 3 /im bump in 
these sources. In this picture, the inner (gaseous) disk is 
by and large optically thin to the stellar radiation, so the 
stellar photons penetrate all the way to Tgub, where they 
are absorbed by the dust. The absorbed radiation heats 
the gas at this location (the rim of the dusty disk) and 
causes the rim to puff up. The details of this model and of 
the refinements that it has undergone since its inception 
are presented in the DM10 review. By increasing the sur- 
face area of the disk material that reprocesses the stellar 
radiation, the puffed-up inner rim naturally contributes 
to the NIR bump emission. However, the estimated ver- 
tical extent of the rim is evidently too small to fully ac- 
count for the strongest ^ 3 / im bumps that have been 
observed in e ither HAe fe.g.. [Tannirkulam et ai]|2008bD 
or cTT (•e.g.. lAkeson et al.ll2005D stars. 

The advances made in infrared interferometry over the 
last decade have led to a powerful new tool for prob- 
ing the origin of the NIR excess in protostars. One 
important early result has been the inference that the 
radius of an "emission ring" model for the interferomet- 
ric visibility data scales as the square root of the stellar 
luminosity i» over four decades in luminosity covering 
cTT and HAe stars as well as l ow-luminosity HBe stars 
(jMonnier fc Millan-Gabetl[200l see Figure 7 in DM10). 

^ ITambovtseva fc Grininl I I2008I) also demonstrated that a dusty, 
centrifugally driven disk outflow can intercept and reprocess a sig- 
nificant fraction of the stellar radiation in this class of sources, but 
they did not explicitly calculate the resulting infrared spectrum. 



This "size-luminosity" relation is consistent with the ex- 
pected scaling of Vguh with (see Equation (|6])), and 
the normalization of the measured correlation in fact im- 
plies a temperature that is close to the expected sublima- 
tion temperatures of grains^ This result thus provides 
direct support for the physical picture underlying both 
the "dusty hydromagnetic disk wind" and the "puffed- 
up inner rim" scenarios. However, recent high-resolution 
observations of two bright HAe stars, AB Auriga and 
MWC 275 (= HD 163296), have challenged the basic 
rim picture by revealing the absence of a "bounce" in the 
visibility curves after their initial drop with increasin g 
baseline (iTannirkulam et al.|[2008bt iBenistv et aIll2010D . 
This behavior implies that the NIR emission region in 
the corresponding sources is not confined to a narrow 
rim but, rather, is spatially extended. 

In this paper we adopt the dusty disk-wind interpreta- 
tion of the NIR excess emission in protostars and develop 
detailed diagnostic tools that enable us to apply it to ob- 
servational data. In Section [2] we review the physics of 
centrifugally driven disk winds and describe our model- 
ing setup. We also briefly describe the Monte Carlo ra- 
diative transfer (MCRT) code that we constructed for 
this project, for which we provide full details in Ap- 
pendix [BJ In Section [3] we present our model results for 
the spectral energy distribution (SED) and the visibility 
curve for several representative wind models and mass 
outflow rates. We compare the results with the data for 
the two aforementioned HAe stars, which have been ob- 
serve d with an unprecedented s ub-milliarcsecond resolu- 
tion ([Tannirkulam et al]|2008a|) . We demonstrate that, 
unlike the "puffed-up inner rim" model, the hydromag- 
netic disk-wind model can in principle account for both 
the magnitude of the ^ 3 /im bump and the shape of the 
visibility curve in these sources without having to invoke 
any other emission component. We discuss additional as- 
pects of this problem in Section [31 where we also further 
comment on the possible relationship between protostel- 
lar and AGN dusty disk outflows, and we then summa- 
rize in Section [5] The extension of this work to higher- 
luminosity protostars, for which radiation-pressure ef- 
fects on dust must also be taken into account, will be 
presented in a future publication. 

2. FORMULATION 

2.1. Centrifugally Driven Winds 

Hydromagnetic driving is thought to be the most likely 
mechanism of accelerating the powerful outflows ob- 
served in protostars and AGNs in view of the fact that 
alternative mechanisms (in particular, thermal and ra- 
diative driving) are generally too weak to account for 
the large momentum and energy discharges inferred in 
these sources. This mechanism is based on the presence 
of a large-scale, ordered magnetic field that threads the 
accretion disk that surrounds the central mass. In the 
context of protostars, this field can be naturally iden- 
tified with the interstellar field that permeates the na- 
tal molecular cloud core and that is dragged in by the 

^ The scaling rc^uh reflects the inverse-square depen- 

dence of the radiative flux on distance from the source and arises 
from the fact that the location of the sublimation radius for given 
dust properties is de termined solely by the incident flux (e.g., 
Ilvezic fc Elitzuri[T997fl . 
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infalling gas once gravitational collapse sets in. Alter- 
native possibilities are that the field is generated by a 
disk dynamo or, as noted in Section [1] that it originates 
in the star. Magnetic stresses could launch outflows in 
different ways, but the most commonly invoked process 
is centrifugal driving, wherein disk material is flung out 
along magnetic field lines that are inclined at a suffi- 
ciently large angle to the disk surface — in analogy with 
the motion of beads along rotating, tilted, rigid wires (see 
iBlandford &: Pavnelll982( ). As the outflowing gas climbs 
above the disk surface, the gradient in the magnetic pres- 
sure associated with the azimuthal field component fur- 
ther accelerates the flow. The disk wind is collimated 
by the hoop stress exerted by the azimuthal field com- 
ponent as well as by the magnetic tension force acting 
along the poloidal field component. Detailed accounts of 
the propert ies of such outflows can be found in the liter- 
ature (e.g., IKonigl fc Pudritzl [20001: fPudritz et all [20071 : 
IKonid fc Salmeronll201lD " 

2.2. Model Setup 

All in all, centrifugally driven winds can be acceler- 
ated and collimated very efficiently, leading to a den- 
sity distribution above the disk t hat is strong ly strati- 
fied in the vertical direction (e.g.. ISafi"eillf QQSal Tbl). Such 
winds are generally also efficient at transporting angular 
momentum away from the disk surface. In fact, there 
is a growing number of protostellar outflows in which 
measured transverse velocity gradients, interpreted in 
terms of rotation in a centrifugally driven wind, indi- 
cate a disk wind that originates at a radius > 1 AU 
and transports a signiflcant fraction of the disk angu- 
lar momentum in the la unching reg ion fe.g.. [ Ray et al.l 



[20071: iChrvsostomou et a l. 2008: Co ffey et al.ll2008[ ). The 
inferred values of the wind launching radius are consis- 
tent with the notion that the outflows can extend beyond 
the dust sublimation radius, which for bright HAe stars 
and their inferred dust properties (grain size 1 ^m; 
sublimation temperature ^ 1850 K; see Section [3]) has 
a magnitude rguh ~ 0.2 (i,/40 L©)^/^ AU (see Equa- 
tion (|6])). There is still no direct evidence for the uplift- 
ing of dust from the inner regions of protostellar disks, 
but se veral different obs e rvations (e.g.. iNeckel fc Staude 
GueFh et al.|[200a iNisini et al1 l2005t iSmith et al. 
Chrvsostomou et al.l 120071 ) provide indirect sup- 
port for this picture. In the case of cTT and HAe 
stars, there are also indications from hi gh- resolution 
spect r o-interferometric o b servations (e . g.. iTatulli et al.l 



1995; 



2005 



[20071 iKraus et al.l [200l [Eisner et ahl |2010D for B^- 

cmitting outflows that originate within Tgub- 

In concordance with these flndings, we adopt a model 
in which a centrifugally driven wind is launched between 
an inner radius, 7'omin, whose value is set to be 0.05 AU 
in all of our models, and an outer radius, romax, which 
we take to lie at 2 rgub (where rgub depends on the source 
parameters; see Equation ((6))). A schematic representa- 
tion of this model is shown in Figured] The inner radius 
could (but need not) correspond to the disk truncation 
radius, which, in turn, could be determined by the stel- 
lar magnetic field. There have already been several re- 
ported measurements o f a kG-strength dipolar magnetic 
field in HAe stars (e.g.. lAlecian et al.|[2009l: iHubrig et al.l 
[201 If) , although it is still unclear whether such a rela- 



tively strong field is present in only a small fraction of 
these sources (e.g.. lWade et al.ll201lD . Using typical val- 
ues for the mass accretion rate (Min) and the stellar mass 
(Af*) and radius (i?*), the steady-state magnetospheric 
truncation radius corresponding to a dipolar field of 
equatorial surface strength can be estimated to be 
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(e.g., [Long et all I2005[) . Our choice of rpmax is mo- 
tivated by the arg ument (e.g., ICasse fc Ferreiral [2000l: 
iPesenti et al.l[200^ that the temperature near the disk 
surface has to be high enough to ensure adequate mass 
loading of the wind. One likely source of heating for the 
gas in this region is stellar radiation, particularly FUV 
and X-ray photons, which (assuming a standard ISM gas- 
to-dust ratio) will reach the disk surface provided that 
the intervening hydrogen-nu cleus column density doe s 
not exceed 10^^ cm"^ (e.g.. lGorti fc Hollenbach|[2009D . 
As we demonstrate in Appendix [XJ stellar photons that 
reach the disk surface penetrate (for the typical parame- 
ters of our disk-wind models) through a dusty wind col- 
umn > 10^^ cm~^ over a distance (along the disk surface) 
^ ^sub, which suggests that a robust outflow is unlikely 
to persist beyond ~ 2rsub- 

To carry out the radiative transfer calculation in the 
dusty wind, we need the gas density structure of the 
outflow and the dust distribution within the wind. For 
the density st ructure , we us e the analytic approxima- 
tions given in iSafied (|1993bl: see his Appendix A) for 
several representat ive wind solut ions that were derived 
semianalytically in Safied (ll993aD using the approach of 
IBlandford fc Pavnd (119821 ). These solutions assume a 



steady, axisymmetric, effectively cold, and radially self- 
similar disk outfiow. The self-similar wind model strictly 
applies to an infinitely extended disk, but, given that the 
density structure changes mainly in the vertical direction 
on a scale that is much smaller than the spherical radius 
i?, and that our results are not sensitive to the values 
of romin and romax, we feel justified in employing them 
in our calculations. For the dust distribution, we simply 
assume a constant gas-to-dust mass ratio of 100 through- 
out the outflow. Although this ratio is expected to vary 
with location in a real disk wind, a more elaborate for- 
mulation is not warranted here in view of the fact that 
our spectral calculation also employs the approximation 
of single-size grains (see Section 12. 3p and that grains of 
one particular size ('~ 1 /.tm) indeed provide a good flt to 
the measured SEDs (see Section [XT]) [1 

One reason for why a constant gas-to-dust ratio is unlikely to 
apply in a real disk outflow is that, for given wind parameters, 
there is an upper bound on the size of grains that can be uplifted 
from the disk at any given radius (see Equation lO). Further- 
more, even grains that are dragged above the disk surface may 
stall and eventually fall b ack in if they are not small enough (e.g., 
ISalmeron fc Irelandll20l3) . In addition, the fact that larger grains 
are characterized by a higher sublimation temperature (see Equa- 
tion ^ and Table 2) implies that the surface layers of the disk in 
the vicinity of r^uj, are preferentially populated by comparatively 
large grains. More detailed future calculations should bo able to 
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Sublimation Radius r^ut, Dust-Free Disl< 

Fig. 1. — Schematic representation of a centrifugally driven disk wind from the vicinity of a low-luminosity protostar. The gaseous wind 
uplifts dust from the disl? surface beyond the sublimation radius rgui, but remains dust-free closer to the protostar. The innermost radius 
^'Omin of the wind-launching zone may (but need not) coincide with the radius rm (Equation |(TJ) where the disk is truncated by the stellar 
magnetosphere. The outer radius romax of this zone is possibly limited to ~ 2 r^ub by dust absorption of stellar FUV photons in the wind. 
Most of the NIR excess emission is produced by reprocessing of the stellar radiation by the dusty wind. 



In the radial self-similarity formulation, any point in 
the wind is specified by the radius vector R, whose cylin- 
drical coordinates {r, (f> z} are 

R=Ne(x),0, rox}- (2) 

Here x = z/ro is the dimensionless vertical coordinate of 
a point along a streamline whose footprint intersects the 
disk at a radius ro (with the subscript '0' denoting the 
equatorial plane). All physical quantities can be repre- 
sented as products of a power law in tq and a function 
of X. In particular, the mass density can be written as 



P = Pi 



-3/2 



vix) 



(3) 



where pi is the density at the disk surface at the fidu- 
cial disk radius ri , which we henceforth take to be 1 AU, 
and ri{x) is obtained from the solution of the MHD wind 
equations. Using Eciuatio ns (2.1), (2.7c), and (5.2) in 
iBlandford fc Pavnd (jl982| ). and normalizing by values 
appropriate to our modeled HAe stars, we can express 
the fiducial density as 



Pi 



=3.5 X 10 



-15 



Mo 



I 5 X 10"** Mq yr- 

IpO In (romax/f Oinin) ' 



2 Ma 



-1/2 



0.03 



2.5 



g cm 



(4) 



where Mout is the total mass outfiow rate from the disk 
(between rgmin and romax) and ipo is the ratio of the verti- 
cal speed Vz to the Keplerian speed v k = {G M^JroY^"^ at 
the disk surface (cf. Equation (20) in lSafi er 1993;a!). The 
radially self-similar wind solutions are defined by three 
parameters: k, the normalized mass-to-magnetic flux ra- 
tio, which characterizes the mass loading of the wind; A, 
the normalized total (particle and magnetic) specific an- 
gular momentum, which describes the angular momen- 
tum transport efficiency of the wind; and = tan 
(where is the angle between the poloidal field compo- 
nent and the disk normal at the disk surface) , which mea- 
sures the initial inclination of the magnetic field lines0 

take these effects into account. 

* In order for a cold outflow to be driven centrifugally from the 
surface of a Keplerian disk, must exceed 30° . 



Table 1 lists the parameters of the three representative 
wind solutions from Saficr (1993J3) that we utilize in our 
calculations. 

If the wind transports away all the excess angular mo- 
mentum of the accreted gas, and the mass accretion rate 
through the disk, Min, remains nearly constant with ra- 
dius, then the ratio of M^ut to Mjn c an be approximated 
byhi{(romax/fOmin)/[2(A - 1)]} (see IKonigl fc SalmeronI 
120 111) . According to this estimate, Maut/Mi-a should lie 
in the range ^ 1 — 5 % for the models listed in Table 1. 
This is consistent with the observational determinations 
of this ratio in low-mass protostars. However, using our 
best-fit values of A/out for the two HAe stars that we 
model in this paper (which are > 4 x 10~® Mq yr~^; see 
Figure the above expression would imply mass ac- 
cretion rates that are in excess of those that have been 
inferred observationally for these two sources (which are 
< lO~^M0yr~^; see below). The discrepancy could be 

reduced if we have overestimated A/out (for example, be- 
cause of our neglect of other likely contributions to the 
NIR excess emission, as discussed in Section[31 or because 
we have overestimated the ratio rgmax/^'Omin in Equa- 
tion(|3])). Note in this connection that a mass outfiow 
rate ~ 1 x 10"® Mq yr"^ was inferred in MWC 275 from 
observations of Herbig-Haro knots in the bipolar outfiow 
emanating from this source (jWassell et al.l [20061 ) . It is 
also conceivable that the mass accretion rate in these 
intermediate-mass sources undergoes a measurable de- 
creases with radius because of a considerable mass loss 
to the wind, so that M\-a_, which is determined from ob- 
servations near the stellar surface, could be significantly 
smaller than the local mass accretion rate near rgub- 
(Thi s situation is thought to po ssibly occur also in AGNs; 
e.g., lElitzur fc Shlosmanll2006l ) Intrinsic variability in 
the source could further contribute to the apparent dis- 
crepancy. For example, MWC 275 exhibited a NIR flux 
increase that did not simply correspond to the behavior 
of accretion diagn ostics such as the Ca II triplet lines 
(jSitko et al.ll2008D . If the NIR variability is associated 
with changes in the outflow rate near rsnh, which could, 
in turn, be related to localized variations in the accretion 
rate at that radius (see SectionH]) , then it is entirely plau- 
sible that the inferred values of Mont and of Mi^ (which, 
as we noted above, is determined at a much smaller ra- 
dius) would not obey the above scaling relation even if 
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TABLE 1 
Wind Solution Parameters 



Model 


K 


A 


5.', 


Notes 


C 


0.01 


75.43 


1.73 


relatively fast wind, compar- 
atively inefficient mass load- 
ing, narrow launching angle 


E 


0.10 


25.63 


3.73 


relatively slow wind, efficient 
mass loading, wide launching 
angle 


G 


0.01 


189.34 


3.73 


fastest wind model, compar- 
atively inefficient mass load- 
ing, wide launching angle 



they did in a strictly steady state. 

Apart from the preceding considerations, the deter- 
mination of the mass accretion rates in HAe and HBe 
stars is complicated by the fact that these objects are 
intrinsically bright at the same UV wavelengths where 
the emission from an accretion shock on the stellar sur- 
face is expected to peak. Furthermore, in contrast with 
cTT stars, it is not clear that the commonly invoked 
magnetospheric accretion model is generally applicable 
to the higher-mass protostars, and there is also a rela- 
tive dearth of metal absorption lines available for mea- 
surements of veiling (an accretion diagnostic). In addi- 
tion, Br7 emission, which is often used to infer Mjn in 
lower-luminosity protostars, can have a significant contri- 
bution from an outflow component in many HAcB e stars 
(e.g. [Tatuni ct al., 2007; Kraus ct al. 2008. ; Wcigch^ et al.l 
120111) . The prevailing uncertainties are reflected in the 

range of values of M-m that have been estimated for 
the two HAe stars under consideration. For exam ple, 
in the case of MWC 275, iGarcia Lopez et al.l ()2006D in- 
ferred a value of 7.6 x 10~* Mq yr~^ based on a mea- 
sureme nt of the Br7 line luminosity. However, even 
though lEisner et al.l ()2010[ ) measured a similar Br7 lu- 
minosity, they deduced a higher value of (~ 1.75 x 
10~^ Mq yi~^) on account of their use of a d ifferent inter- 
pretive scheme. iDonehew fc BrittainI ()2011|) . in turn, in- 
ferred Min ~ 6.9 X 10~^ Mq yr~^ for this source by using 
the so-called Balmcr discontinuity as a veiling tracer, but 
[Mcndigutfa ct al . (2011,), employing a similar approach, 
deduced only an up per limit of ^ 3.1 x 10~^ M r:^ yr~^ . In 
the case of AB Aur. IGarcia Lopez et a l.' (200&) estimated 
Min ^ 1-4 X 10""^ Mq yr~^, whereas iDonchcw fc BrittainI 
(|2011D inferred a value of only - 1.8 x 10"^ A/oyr-^0 

2.3. Spectral Calculations 

All the spectral results presented in this paper were 
obtained with the help of a newly developed and fully 
tested MCRT code that is based on the radiative equi- 
librium and temperature correction scheme discussed in 
iBiorkman fc Woo"dl ()2001h : it is described in detail in Ap- 
pendix [BJ Although this code is technically 2D, it fully 
captures the actual 3D trajectories of the modeled pho- 
tons under the postulated axisymmetry. In any given run 
of the code, the grains are assumed to have one particu- 
lar size: either small (0.001 /xm), i nterme d iate (0 .1 /xm), 
or large (1.0 ^m). As discussed in iSafieii (|1993a[) . there 

^ Note in this connection that, by attributi ng the origin of the 
hydrogen infrared lines in AB Aur to a wind, INisini et al.l 1)19951 ) 
estimated a mass outflow rate of ~ 3.3 x lO^^MQyr"-*^ in this 
source. 



is a maximum size Cmax of grains at any given radial lo- 
cation for which the wind has sufficient momentum to 
overcome tidal gravity and uplift the dust from the disk 
surface. Combining Equation (E7) in that paper with 
Equation (|4]) above, we estimate 

''--^''^ " (sTTF^So^) (^) (^) 

2.5 ) VO.SAUJ ' 

(5) 

which indicates that grains much larger than ~ 1 fim are 
unlikely to be present in typical disk outflows from HAe 
sources. (Although our model neglects dynamical inter- 
actions of dust particles in the wind, such interactions 
are not expected to lead to measurable growth after the 
grains enter the outflow.) 

We assume that the dust consists of pure silicates and 
we employ grain optical properties (derived under the 
assumption that the grains are compact and spherical) 
from the tabulation of iWeingartner fc Draini (|200lD FI 
These properties depend on the grain size. Under con- 
ditions of radiative equilibrium, a given grain emits as 
much energy per unit time as it absorbs. The dust 
optical properties at the wavelengths where the stellar 
emission peaks determine how effectively this radiation 
is absorbed, whereas the corresponding properties at the 
longer wavelengths where most of the dust emission oc- 
curs determine how effectively the absorbed radiation is 
reemitted. Grains of size a <^ 1 /im typically have a low 
cooling efficiency, whereas grains whose size exceeds a 
few microns have an effectively gray opacity and reradi- 
ate the absorbed stellar radiation with a comparatively 
high efficiency (see DM10). Since the radiative equilib- 
rium condition determines the dust temperature, large 
grains tend to be cooler than small grains for a given 
incident flux. By the same token, large grains attain a 
given equilibrium temperature at a smaller radius (where 
the incident stellar flux is higher) than small grains. This 
is quantifled by the scaling function H in the expression 

r,nb = i/(r*,r,nb,a)(L*/L0)i/2r-2 AU (6) 

for the sublimation radius, which is obtained from the 
radiative equilibrium condition by substituting Tgub for 
the dust temperature. Table 2 lists the values of 
H{T^,Tsuh, o-) for our representative grain sizes and for 

® Available at www.astro.princeton.edu/~draine/dust/ 
dust.dicl.html. 
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three values of the dust subhmation temperature, assum- 
ing that the stehar photosphere is characterized by an 
effective temperature T, « 9650 K. 

The dust subhmation temperature depends on the 
grain composition as wcU as on the ambient gas pres- 
sure. Recent apphcations in the hterat ure have often 
used t he empirical relation presented by llsella fc Nattal 
(l2005l) . which gives Tsub for silicates as a function of the 
lo cal gas density. Th is relation is based on the results 
of lPoUack et al.l (|1994l ). who fitted a Clausius-Clapeyron 
equation to data for Olivine. The extraction of a den- 
sity dependence from this equation was done under the 
implicit assumption that the dust and gas temperatures 
are equal, but this is typically not the case for the com- 
paratively low-density disk winds that we consider. In 
view of this, and in the interest of simplicity, we neglect 
any density dependence of the sublimation temperature 
in our treatment and content ourselves with using a sin- 
gle sublimation temperature (having one of three rep- 
resentative values) in any given runQ Given that pre- 
vious inferences of the sublimation temperature, made 
in the context of the in ner rim model, have ranged from 
- 1150 K for MWC 275 (iBenistv et aLllMol) to - 1950 K 
for AB Aur (jTannirkulam et al.l l2008b[) . we adopt the 
following representative values for Tgub^ 1250 K, 1500 K, 
and 1850 K. As Table 2 demonstrates, the parameter de- 
pendence of Tsub is not, in general, monotonic: For the 
given value of T, and with Tguh fixed, the value of Tgub 
peaks for the intermediate-size grains. 

In our formulation, we represent the star in a simpli- 
fied fashion as a point source at the origin of the coordi- 
nate system. We further simplify the treatment by tak- 
ing the disk to be a flat surface in the equatorial plane. 
In this limit, no stellar photons impinge on the disk di- 
rectly. Radiation does, however, reach the disk as a re- 
sult of scattering and reemission of stellar photons in the 
dusty wind. We treat the interaction of this radiation 
with the disk in a manne r analogous to the two-zone ap- 
proximation described in |Chi^g[^^oIdreic3 (jl993). In 
this picture, the incident radiation is fully absorbed in 
a dusty surface layer whose temperature is calculated in 
the same way as that of the grains in our wind model. 
The heated surface layer reemits the incident radiation 
at longer wavelengths, and it is assumed that half of this 
radiation is directed outward and the other half is ab- 
sorbed and thermalized in the disk interior, which then 
reemits it as a blackbody. We neglect, however, any po- 
tential contribution to the disk emission from energy dis- 
sipated in the accretion process. Our MCRT code follows 
the paths of photons that are emitted by the disk in the 
same manner that it handles the transfer of the stellar 
radiation through the wind. 

In this paper we concentrate on the dusty wind and 
disk region beyond the sublimation radius, and our spec- 
tral calculations initially neglect any contribution of the 
disk or the wind at smaller radii. We nevertheless return 
to address this issue in Section [3.31 where we estimate in 
an approximate manner the effects of the thermal emis- 
sion from the interior wind (due mainly to intrinsic am- 

We have, however, verified that our best-fit models do not sig- 
nificantly change if we employ instead a density-d ependent subli- 
matio n temperature using the expression derived bv llsella fc Natt"al 



TABLE 2 
Scaling function H in Equation (O 





a 


(^tm) 






0.001 


0.1 


1.0 


1250 


2.18e5 


2.95e5 


1.23e5 


1500 


2.15e5 


2.85e5 


1.15e5 


1850 


2.05e5 


2.65e5 


1.09e5 



bipolar diffusion heating) and disk (due to reradiation of 
directly impinging stellar photons). We do not, however, 
include the effect of "back warming" of the interior disk 
and wind region (or of the star) by photons that originate 
beyond r^nh- 

3. RESULTS 

In this Section we present the results of our spectral 
calculations, employing different disk-wind models, mass 
outflow rates, grain sizes, and dust sublimation temper- 
atures to infer the dependence of the NIR observational 
signatures on the relevant physical parameters and to de- 
duce the parameter values that best fit the data for AB 
Aur and MWC 275. We first consider the dusty wind 
and disk region at r > rgub and describe, in turn, the 
predicted spectral energy distributions (SEDs) and in- 
terferometric visibility curves , and we then discuss the 
expected contribution from the gas in the disk and the 
outflow interior to Vsuh- 

3.1. Model SEDs 

As was alread y found in previo us studies (e.g., 
llsella et al.l l2006t iTannirkulam et al.l l2008b[ ). the dust 
properties that most directly affect the appearance of the 
NIR bump are the size of the grains (a) and the sublima- 
tion temperature (Tgub)- Besides these two parameters, 
we consider also the dependence of the calculated SEDs 
on the specific disk-wind model (C, E, or G; see Table 1) 
and on the mass outflow rate (Afout)- To fully specify a 
solution, we also need the values of the stellar mass (Af*) 
, total luminosity (i*), and effective temperature (r*), 
which are inferred observationally for any given sourceQ 
as well as the values of the parameter ipQ and of the ratio 
fomax/^'omin (scc Equation ([4])) — we flx the value of ipo 
to be 0.03 and evaluate the ratio of the outermost and 
innermost wind radii using the prescription outlined in 
Section [2^ 

Figures [2] and [3] show results from our model calcula- 
tions as well as data points for AB Aur and MWC 275. 
Since our interest centers on the NIR spectral regime, 
the SEDs arc plotted only up to wavelengths of ^ 10 fim. 
It is immediately apparent that this model can account 
for the prominent NIR bumps exhibited by these sources 
while also roughly reproducing the overall shapes of the 
measured spectra. The main contribution to the NIR 
excess emission comes from the dusty wind, which in- 
tercepts the stellar radiation in a "reprocessing surface" 
whose effective extent (for a given wind model) depends 
mainly on the mass outflow rate. While this surface is 

° We adopt the values given in lGarcia? Lopez et aLl 120011 : M* = 
2.4 Mq, T. = 9840 K, L, = 48 Lq for AB Aur, and M, = 2.3 Mq, 
T* = 9450 K, L* = 36 L0 for MWC 275. 
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Fig. 2. — SED fits for AB Auriga, showing the effects of varying different model parameters. In each of the four panels, the solid (red) 
curve represents the base model, whereas the dashed (green) and dotted (blue) curves show the changes to the predicted spectrum that arc 
induced by varying the value of one of the parameters: the wind model (a), the grain size (b), the sublimation temperature (c), and the 
mass outflo w rate (d). The d ata points are from the MDM telescope (filled circles; Tamiirkulam et al. 2008b) and from ISO [thick solid 
(gray) line: IMeeus et aT1l2001l and references therein]. The stellar emission is represented by a blackbody spectrum (thin dotted line). 




Fig. 3. — Best-fit curves to the SEDs of AB Auriga (a) and MWC 275 (b). In the case of AB Aur, the two curves demonstrate the 
possible parameter degeneracy that can be encountered by fitting the SED data alone. In the case of MWC 275, the two curves (which 
both correspond to Tguj, = 1500 K) demonstrate how a change in the wind outflow rate could potentially accoun t for the NIR vari ability 
exhibited by this source. The high-state (squares) and low-state (triangles) measurements for MWC 275 are from lSitko et all l|2008I V The 
references for the data points from the MDM telescope and from ISO are given in the caption to Figure [2] 



analogous to the puffed-up inner rim, it can be signif- 
icantly larger for a sufficiently high value of Mout; the 
wind model can therefore explain even the strong ~ 3 /xm 
bumps seen in these bright HAe sources, for which the 
rim model evidently falls short. Our spectral model also 
includes a dusty disk component, which begins to make a 
significant contribution to the spectrum at wavelengths 
> 4 /xm. Since our disk model is rather simplified (es- 
pecially in its assumption of geometrical flatness), the 
predictions at those longer wavelength are not expected 
to be accurate. As the spectral regime beyond the NIR 
bump is not the focus of this work, we have not at- 
tempted to improve on our disk emission model, but we 



note that more realistic models of passive flared disks 
can r eproduce the observ ed spectra at these wavelengths 
(e.g., IChiang et al][200l . 

The dependence of the synthetic spectra on the differ- 
ent model parameters is illustrated in Figure [2] with fits 
to the data from AB Aur. The base model corresponds 
to 1 fiia grains with a sublimation temperature of 1500 K 
that are embedded in a wind that has a mass outflow rate 
of 7 X 10~^ Mq yr~^ and is described by wind model C, 
and each of the four panels presents variations in one 
of these attributes. Panel (a) shows the dependence of 
the spectra on the choice of wind model. The domi- 
nant effect on the emitted flux level appears to come 
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from the self-similarity parameter ^q, which character- 
izes the initial opening angle of the poloidal streamlines 
(see Section 12. 2p : The NIR flux is measurably higher 
in model C, which has the narrowest launching angle, 
but there is not much difference in the spectral prop- 
erties of models E and G, which have disparate values 
of the self-similarity parameters k and A but the same 
(higher) value of (see Table 1). Panel (b) examines 
the dependence on the grain size a. As was noted in 
Section 12. 3[ large {a > 1 fim) grains cool with a com- 
paratively high efficiency, which explains why the wind 
contribution to the ~ 3 fiia bump increases with a. The 
large grains, with their effectively gray opacity, also bet- 
ter match the flat shape of the spectrum in the bump 
region. It was further noted in Section [2731 that the value 
of the sublimation radius is largest for the intermediate- 
size (a ~ 0.1 /im) grains (sec Equation ([6|) and Table 2). 
This implies that, for these grains, the emission at any 
given wavelength originates further out in the disk than 
it does for the other grain sizes. At those larger distances, 
the wind density is lower because of the r~^/^ scaling of p 
(sec Equation ([3])) and of the inverse linear dependence 
of the normalization pi on ln(romax/''Omin) (see Equa- 
tion (j4|)). Therefore, the fraction of the stellar radiation 
that is intercepted by the wind and scattered down into 
the disk is lower in this case, resulting in a cooler disk 
that contributes more to the flux at longer wavelengths. 
This explains why the intermediate-size grains exhibit 
the strongest ~ 10 pm emission feature. 

Panel (c) of Figure [2] depicts how the predicted spectra 
are modified when the value of Tgub is changed. Since a 
higher sublimation temperature implies a smaller value of 
the sublimation radius, the dusty emission region moves 
inward to where the wind density is higher (see equa- 
tions ([3]) and ([4|)), leading to a higher optical depth in 
absorption and hence to more reprocessing of the stellar 
radiation. Therefore, we can expect stronger NIR emis- 
sion for higher values of Tgub- This expectation is borne 
out in the figure, which shows that the flux becomes pro- 
gressively higher as Tguh is increased. A similar trend is 
expected when the mass outflow rate goes up, in view of 
the linear dependence of the density on Mout (see Equa- 
tion ([3])). This expectation, too, is confirmed by the 
explicit calculation, as is shown in panel (d). 

Figure [3] presents our "best-fit" SED models for AB 
Aur and MWC 275. Because of the limited number of 
runs of the MCRT code that we have performed, our 
chosen parameter values are only accurate, on average, 
to ~ 10 %. Furthermore, as has already been known from 
previous work in this area (see DM10), a fit to the SED 
docs not uniquely fix the underlying physical model. We 
illustrate this fact within the framework of our scheme in 
panel (a) of this figure, where we show two comparable 
fits to the spectrum of AB Aur that are based on differ- 
ent combinations of the parameters Mout and Tguh- As 
we demonstrate in Section 13. 2| this parameter degener- 
acy can be lifted by considering also the best fit to the 
visibility data. 

Panel (b) of Figure [3] shows SED flts for MWC 275. 
This source exhibited a ^ 30 % variabihty in the 1 — 
5^m flux level around 2002 (jSitko et al.ll2008| ). and the 
data points associated with this "outburst" (the "high" 
state) and with the quiescent emission (the "low" state) 



are indicated in the flgure. We fit the "high" state of 
this source with the base model of Figure [U whereas an 
approximate match to the 'low" state can be obtained by 
reducing the mass outflow rate in this model by Ail/out ~ 
3 X 10~^ Mq yr~^. Physically, the outburst exhibited by 
this source could be associated with the injection of a 
density inhomogeneity (a "cloud" or "clouds") into the 
wind. Such clouds would move across the line of sight at 
approximately the wind speed and could in principle also 
account, through shadowing, for the variability that was 
measured in this source in the scattered-light image of 
the outer disk (whose spatial sc ale is a few hundred AU ) 
on a time scale of a few years (jWisniewski et al.l[200l . 
Further discussion of this scenario is given in Section |4l 
Our model flts in this subsection are based exclusively 
on the contributions of the wind and disk beyond rgub- 
If the inner disk is optically thin to the stellar radiation 
then an additional contribution from the heated inner 
edge of the dusty disk, as envisioned in the "puffed-up 
inner rim" scenario, could be expected. The contribu- 
tion of the latter component to the NIR excess in these 
sources would be subordinate to that of the dusty wind, 
but it could conceivably help improve the detailed flt to 
the measured spectra. A further contribution to the ob- 
served spectrum could come from the disk and wind at 
f < ''sub, as discussed in Section [3l3l Any additional con- 
tribution of this type to the NIR excess emission would 
reduce the inferred value of A/out for the wind compo- 
nent. 

3.2. Synthetic Visibilities 

As was already remarked in Section the ability to 
flt the SED is not sufficient for distinguishing between 
competing interpretations of the NIR excess. In fact , 
spherical dusty halo models (e.g., iVinkovic et al.|[2006[) 
can also account for the NIR bump and, as pointed out 
by DM10, may even be more consistent than the "puffed- 
up inner rim" model with the apparent lack of a clear 
correlation between the NIR flux and the disk inclina- 
tion among the observed sources0 However, with the 
recent advent of NIR intcrferomctry it has become pos- 
sible to start spatially resolving the NIR emission regions 
of HAe and HBe stars and thereby to better discriminate 
among the different models. As was noted in Section [1] 
high-resolution measurements of this type are already 
available for the two HAe stars mo deled in this paper. 
The existing data for these sources (jTannirkulam et al.l 
HoOSb; Bcnisty et al. 20l3) indicate both that a spheri- 
cal halo model is very unlikely to reproduce the observa- 
tions and that (as explained below) a "puffed-up inner 
rim" model is ruled out (unless an additional emission 
component is present within the sublimation radius). In 
this subsection we demonstrate that, in contrast, the in- 
terferomctric results are consistent with the disk-wind 
model. 

While direct imaging is still not technically feasible 
for the inner regions of protostellar systems, the NIR 
emission zones of HAe stars are already being probed 

^ Although the ccntrifugally driven wind model gives rise to a 
nonspherical density distribution, the degree of flattening of this 
distribution is much smaller than that of a thin accretion disk and 
therefore a strong correlation between the reprocessed flux and the 
inclination of the symmetry axis is also not expected in this case. 
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Fig. 4. — Model NIR visibilities for AB Auriga (a) and MWC 27 5 (b), corresponding to the "best fit" parameter combinations presented 
in Figur es [3^ and[3l3, respectively. The VLTI data points arc from lBenistv et al.l 120101 ). whereas the references for the other data sets are 
given in ITannirkulam et al.l Il2008bl) . 



by long-baseline interferometry. In this technique, the 
light from two (or more) telescopes is combined to gen- 
erate an interference patten ("fringes"). The amplitude 
and phase of these fringes, which depend on the pro- 
jected distance on the plane of the sky between any given 
pair of telescopes (the "baseline"), are combined into a 
complex quantity, V (the "visibility" ) . According to the 
van Cittcrt-Zernikc theorem, the visibility represents a 
Fourier transform of the spatial intensity distribution, 

y(w,w) = J J I{a,/3)e^'''^°"'+^'"Uad/3, (7) 

where /(a, (3) is the intensity distribution on the sky as 
a function of the angular coordinates (direction cosines) 
a and (3. The spatial frequencies u and v are the nor- 
malized (by the observation wavelength A) projections of 
the baseline on a plane perpendicular to the line of sight: 
u = Bx/\ and v = By/X. It should be evident from 
the form of Equation ([7| that an extended object with a 
Gaussian-like intensity profile would produce a visibility 
that is also Gaussian-like — smooth and single-peaked. 
By contrast, the Fourier transform of the intensity profile 
of a uniform-brightness ring exhibits the opposite behav- 
ior, displaying large secondary peaks. Thus, the presence 
or absence of a secondary bounce in the interferometric 
visibility curve indicates whether the emission region is 
"sharp" (like the uniform ring) or "fuzzy" (like a Gaus- 
sian) . For a review of v isibilities and of how to interpret 
them, see IBerged (l2003t see also Figure 5 in DM10). 

The basic rim model, by construction, involves a very 
"sharp" emission region and therefore implies a pro- 
nounced secondary bounce in the visibility curve. As we 
already noted in Sectionfl] such a bounce is clearly absent 
in the measured visibilities of AB Aur and MWC 275, 
which are found instead to be relatively flat beyond their 
initial dropoff. We now show that the smooth, extended 
dust distribution that characterizes the hydromagnetic 
disk-wind model naturally accounts for the shape of the 
visibilities measured in these sources. The predicted vis- 
ibilities are derived by carrying out the integral in Equa- 
tion ([T]) using the intensity distribution obtained from 
the MCRT calculation. We consider a wavelength inter- 
val that is centered on 2.2 jim (corresponding to a typical 



K band) and plot all visibilities as a function of the "ef- 
fective" baseline 

Bcs = Bprojoctod [cos^ S + cos^ 6 siu^ J] (8) 

where 6 is the angle between the baseline direction and 
the major axis of the (projected) disk, and e is the angle 
between the dis k's rotation axis and the lin e of sight. As 
pointed out by ITannirkulam et aD (j2008ai rbl) , the quan- 
tity in brackets properly accounts, in the case of a flat 
disk, for the dependence of the interferometric resolution 
on the disk's inclination and position angle. 

Figure |4] shows the visibility curves in the disk-wind 
model that correspond to the "best fit" parameters ob- 
tained from the SED modeling in Section 13.11 Figure 
is for AB Aur and corresponds to the best-fit spectra pre- 
sented in Figure |3]a,, whereas Figure HJd is for MWC 275 
and corresponds to Figure [Jb- In obtaining these curves, 
we used values of the disk inclination angle e of 21° 
and 48° for AB Aur and MWC 275, respectively (see 
ITannirkulam et aD l2008b[ ) . Considering first AB Aur, 
the two curves shown in Figure [3^ correspond to distinct 
combinations of the parameters Afout and Tgub but give, 
as we already noted in Section lO] similarly good fits to 
the SED data for this source. Figure [4^ reveals that the 
visibility data can be used to break this degeneracy: It is 
seen that the Tsub = 1850 K, Afout = 4 x 10"^ A/oyr"! 
model both matches the short- and medium-baseline 
data very well and provides a clearly superior fit to the 

lower- Tsub, higher- A/out model. Figures [3^ and |4^ thus 
illustrate how the simultaneous modeling of SED and 
visibility data can be used to narrow down the choice 
of "best fitting" model parameters. These figures also 
demonstrate that, in contradistinction to the rim model, 
the dusty disk- wind model can by itself account for both 
the SED and the short /medium-baseline visibility data 
for this bright HAe star. However, our model fits under- 
estimate the long-baseline (i?off > 300) visibility data for 
this source. It is quite plausible (in view of the expected 
multiple visibility peaks) that the inclusion of a (subor- 
dinate) inner-rim emission component could improve the 
fit for this spatial-frequency regime, but this remains to 
be demonstrated by an explicit calculation. 
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Unlike the situation depicted in FigureH^, the visibihty 
curves shown in Figure |4}d for MWC 275, which corre- 
spond to the two parameter combinations used in Fig- 
ure [Sb, are barely distinguishable. We recall that these 
two models differ only in the magnitude of Mout, with 
the two listed values chosen to account for the observed 
flux variability in this source. The visibility data shown 
in the figure were collected over different periods and do 
not directly correspond to the SED data presented in 
Figure Eb. The fact that the predicted visibilities are 
nearly identical is consistent with the apparent concor- 
dance of the different data sets in the \ow-BeS regime 
where they overlap. It is conceivable, however, that fu- 
ture multi-baseline, high-time-resolution interferometric 
observations might be able to help identify the physical 
mechanism (e.g., the ejection of dusty clouds; see Sec- 
tion U]) that is responsible for the NIR variability in this 
class of sources. 

Another useful quantity that can be measured interfer- 
ometrically when at least three telescopes are involved 
in the observation is the "closure phase" — the sum 
of the fringe phases for three baselines that form a tri- 
angle. The value of this quantity serves as a probe of 
the degree of point symmetry of the source: The closure 
phase vanishes for a centrosymmetric image, and it in- 
creases the mo re the emission isconcentrated off center 
(or "skewed"). iMonnier et"all (|2006[ ) carried out a NIR 
closure-phase survey of a sample of 14 HAeBe stars (in- 
cluding AB Aur and MWC 275) and found that only six 
(including AB Aur) exhibited statistically significant (al- 
though generally small) values. These results are at odds 
with the original version of the rim model, in which the 
puffcd-up edge of the disk was assumed to be perfectly 
vertical. Since only the side of the rim that faces the 
star emits in this picture, this model predicts NIR images 
that are strongly skewed for nonzero inclination angles. 
The discrepancy is mitigated in th e subsequently devel- 
oped "rounded rirn" mo dels (e.g., Ilsella fc Nattal [20051 : 
iTannirkulam et al.l[2007t see Figure 4 in DM10 for a pic- 
torial representation) . The overall appearance of the NIR 
images that are obtained in the dusty disk-wind model 
closely resembles the basic structure of the images that 
are derived using the rounded-rim model. However, given 
the more extended nature of the emission region in the 
wind model, it may be expected that the NIR images ob- 
tained in this case would be more centrosymmetric than 
those predicted by the rounded-rim model. 

3.3. Contribution of the Region Interior to Vguh 

As discussed in Section [2^ (sec Figure [1}, the wind 
launching region in our model extends down to a radius 
?'Omin that is generally smaller than rgub- As we already 
noted , this picture is supported by Br7 line observations 
(e.g., iTatulli et al] [20071; iKraus et all [2001) . The disk 
truncation radius, which may be determined by the in- 
teraction of the accretion flow with the stellar magneto- 
sphere (see Equation ([T|)), could be even smaller. But 
even though the outflow and the disk at r < Vsuh are 
dust free (unless they contain highly refractory grains), 
they can still contribute to the NIR emission. We now 
estimate this contribution. 

Since the interior outflow does not contain dust, it 
remains optically thin to the bulk of the stellar radia- 



tion. It can, however, be heated by internal processes 
as we ll as by ionizing radiation from the central star. 
iSafieii (Il993al ) examined the thermal structure of disk- 
driven hydromagnetic winds in cTT stars and identified 
a feedback mechanism involving ambipolar diffusion (or 
ion-neutral drag) in the initially weakly ionized disk ma- 
terial that feeds the outfiow. He showed that, for a 
given dynamical structure determined from a solution 
of the wind equations (which is obtained under the usu- 
ally well-justified assumption that thermal forces do not 
strongly infiuence the fiow dynamics), the ambipolar dif- 
fusion heating rate scales inversely with jpi, where 7 is 
the ion-neutral coUisional coupling coefficient and pi is 
the ion mass density. This shows that, s o long a s pi is 
small, the heating could be strong. ISafieii (Il993af ) found 
that this process indeed leads to a rapid increase of the 
temperature to ^ 10'* K (see his Figures 5 and 6): At 
this value the coUisional ionization of hydrogen starts to 
measurably increase pi, which in turn causes the heat- 
ing rate to decline. This mechanism is evidently quite 
robust, although its effectiveness could be reduced if an 
external ionization source is present. We expect the same 
process to operate also in the dust-free regions of HAe 
disk outflows, and in Figure [5] we present a scher natic o f 
the resulting thermal structure based on Safier's ()1993aD 
calculations. The wind can be roughly divided into three 
regions: (1) a cool (~ 1500 K) molecular zone that oc- 
cupies a thin layer at the bottom of the outflow near 
the surface of the disk; (2) a warmer 3000 K) atomic 
and molecular zone that occupies a slightly wider layer 
directly above the cool zone; and (3) a hot (~ 10** K) 
low-density region that lies above the previous two and 
extends to the edge of the modeled domain. In Safler's 
(|1993a| ) calculations, the temperature of the upper layer 
was invariably > 10000 K and it was characterized by a 
high deg ree of ionizat i on. H owever, subsequent calcula- 
tions bv iShang et al.l ()2002D determined that, when the 
temperature is ~ 10000 K and the gas is fully atomic, 
the coupling coefficient 7 is about an o rder o f magn i- 
tudc higher than the "cold" value used in lSafierl ()1993af ). 
which would lead to a commensurate reduction in the 
ambip olar diffusion heating rate. But Shang et al.'s 
()2002|) results also indicate that, so long as the gas is 
mostly molecular, Safier's ()1993a|) adopted value for 7 
should roughly apply (see their Equation (E8)). Given 
that molecular hydrogen is expected to survive in the flow 
until the temperature reaches ~ 5000 K, at which p oint 
it rapidly dissociates (see Section 4.6 in lSafierl[l993a[ ). we 
estimate that the temperature of the uppermost wind re- 
gion would level off at th at value in the abse nce of any 
other heating mechanism. IShang et al.l (l2002f ) suggested 
that one conceivable mechanism for additional heating 
is the dissipation of fluctuations that are induced at the 
source. They wrote down a parametrized phcnomenolog- 
ical expression for the r ate of such mechanical heating, 
aud io 'Brien et all ()2003() showed that, with the inclusion 
of this additional mechanism, a disk wind's temperature 
could in principle again reach (or even surpass) 10000 K. 
For completeness, we thus consider the following two op- 
tions for the temperature T3 of Region 3: T3 = 10000 K, 
for which we label the zone as '3a', and T3 = 5000 K, for 
which we use the label '3b'. Note that, when mechanical 
heating is strong enough to raise the wind temperature to 
10000 K, it may largely eliminate the cooler underlying 
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Fig. 5. — Schematic of the thermal structure of an internally 
heated wind. The two options for the temperature T3 of Region 3 
are discussed in the text. If the wind were heated instead by the 
star's ionizing radiation, its structure would be similar to that of 
Region 3a. 

layers (see Figure 1 in lO'Brien et al.ll2003D . However, as 
we show below, this modification would have little effect 
on the NIR appearance of the source. 
HAeBe stars emit Lyman continuum photons (at a rate 
in the range ~ 10**^ — 10^^ s~^; e.g.. [Alexander et al.l 
|2005( ). which can also contribute to the heating of the 
wind. We can utilize the approximation of a point source 
and th e self-similarity (in the sphe rical radial coordinate) 
of the IBlandford fc Pavni (I1982D wind solution to ex- 
press the outer spherical radius i?out of the ionized zone 
within the wind (the "Stromgren region" ) as a function 
of the polar angle 6*, 

i?max = i?in exp{g, / [47rai?f„7l2 (i?in)] } , (9) 

where Rm (which is also a function of 6) is the spherical 
radius of the intersection point of a ray from the center 
in the direction 6 with the innermost streamline of the 
wind, n(i?in) is the wind's particle density at that loca- 
tion, and a is the hydrogen recombination coefficient; we 
assume that the gas is already atomic in the region of 
interest. Using this expression, we find that only a thin 
layer in the vicinity of the innermost streamline is pho- 
toionized for the values of Mout implied by the spectral 
fitting described in Sections and [X^ We therefore do 
not consider this case further in this Section. We note, 
however, that when Afout is low enough for the ionizing 
photons to penetrate the wind, the shape and temper- 
ature of the Stromgren region closely resemble those of 
Region 3a in Figure [5j 

In calculating the NIR emission from the optically thin, 
dust- free wind component, we used average Planck opac- 
ities that were provided to us by J. Ferguson (private 
communication 2012) since the relevant region in the 
temperature-density parameter space is not included in 
previ ously published opacity tables (e.g.. iFerguson et ahl 
|2005| ). We treated the underlying disk component as a 
flat slab that intercepts the radiation from the star (char- 
acterized by an effective temperature and radius R^) 
and reradiates it as a blackbody of temperature 

/-, \ 1/4 / n \ 1/2 
fe.g- iHubeuA^lflQQOl) . The inner disk should be opti- 



cally thick for the mass accretion rates inferred in our 
model, and w hile this poses a prob lem for the inner rim 
scenario (e.g.. lMuzerolle et al.ll200^ . it has no adverse ef- 
fect on the disk outflow interpretation of the NIR bump, 
in which a dusty wind intercepts the stellar photons well 
above the disk surface. The results wc show were ob- 
tained under the assumption that the disk's emission re- 
gion coincides with the wind launching region (with inner 
radius romin = 0.05 AU). Figure |6] shows the contribu- 
tion of the interior emission components to the SED of 
a source like AB Aur for our two representative choices 
of the temperature in Region 3 of the wind. It is seen 
that the disk emission dominates, but that there is also 
a measurable contribution at the shortest wavelengths 
from wind Region 3, particularly in the T3 = 10000 K 
case. The reason why the main contribution from the 
wind comes from Region 3a is that the emitted flux is 
oc pVT*, and this region has the highest temperature and 
comprises by far the largest volume V, which more than 
make up for its lower density p. For the parameter range 
that we explored, the interior components contribute at 
a level of only ^ 10 — 20 % of the dusty wind's emission. 
However, as Figure |6] demonstrates, this contribution suf- 
fices to reduce the required wind outflow rate by nearly 
a factor of 2 compared to the fit shown in Figure [3^ for 
the same wind model (C) and sublimation temperature 
(1500 K)E3 

Figure [7] shows the effects of the interior emission com- 
ponents depicted in the left panel of Figure [6] on the pre- 
dicted K-band visibility. It is seen that, overall, these 
effects are minor, which is consistent with the relatively 
low contribution of these components to the total flux at 
these wavelengths. The qualitative changes to the vis- 
ibility curve can be understood from the fact that the 
interior components contribute emission at smaller radii: 
This both reduces the effective size of the emission region 
somewhat, which shifts the first minimum to a slightly 
longer baseline, and further "smears out" the emission, 
which is refiected in a small drop in the amplitudes of the 
secondary peaks. The inference that changing the out- 
fiow rate by nearly a factor of 2 in our homogeneous wind 
model does not have a major impact on the visibility 
curve is consistent with the conclusion we have already 
drawn from modeling the visibilities for the "high" and 
"low" states of MWC 275 (sec Figure Sb). The fact that 
the addition of the interior components does not modify 
the predicted visibility in a significant way implies that 
our best-fit model for AB Aur — which is reproduced 
here from Figure|4^ — still requires a comparatively high 
dust sublimation temperaturer^ 

4. DISCUSSION 

The availability of both SED and high-resolution visi- 
bility data for several bright HAe stars has made it pos- 
sible to place strong constraints on physical models of 

We do not use our best-fit model for AB Aur, which is char- 
acterized by Tgub = 1850 K, to illustrate the contribution of the 
interior wind and disk components to the SED because this con- 
tribution becomes hard to discern in the corresponding figure for 
the higher- Tsub case. 

Th is result is similar to that obtained by ITannirkulam et al.l 
Il2008bl) in the context of the rim model. However, in that paper 
they also needed to invoke a major interior emission component to 
account for the full strength of the NIR bump and for the lack of 
a pronounced secondary bounce in the visibility. 



12 



Bans and Konigl 



Inner Wind Region 1 (1500 K) 
Inner Wind Region 2 (3000 K) 
Inner Wind Region 3a (10000 K) 
inner Wind Region 3b (5000 K) 
Interior Disk Emission 



Total Interior Emission 
Dusty Wind Emission 
Dusty Wind + Interior Emission 
MDM Data 
ISO Data 



1e-11 



CM 

E 

~w5 Ie-12 

-I—' 
-I—' 

CO 
LL 

c< 1e-13 



1e-14 



\ • 




_! 1 1 1 I L 



10 



10 



X (|a,m) 



Fig. 6. — Contribution of the interior wind and disk to the model SED fit for AB Auriga. The contribution of the three wind regions 
described in Figure \5\ as well as of the disk component are shown for the two adopted representative temperatures for wind Region 3: 
T3 = 10000 K (Left) and T3 = 5000 K (Right). Also shown are the T^uh = 1500 K, A^out = 4 X 10"* Mq yr"! dusty wind model from panel 
(d) of Figure [5] and the data points displayed in that figure. 
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Fig. 7. — Contribution of the interior wind and disk to the model NIR visibility fit for AB Auriga. The dusty wind model and the interior 
wind and disk components are those shown in the left panel of Figure |6] The data points are the same as those displayed in Figure |4^. 
Also shown for reference is the best-fit model from Figure |4^, corresponding to a sublimation temperature of 1850 K and an outflow rate 
of 4 X lO-^Mgyr-i. 
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the NIR excess emission in these sources. The observa- 
tions are consistent with reprocessing of the stellar ra- 
diation by a nonspherical distribution of dust, and they 
have been commonly interpreted in terms of a disk with 
a puffed- up inner rim. The flat shape of the spectrum in 
the ^ 3 fim bump regime points to comparatively large 
(> 1 /im) grains that have a near-gray opacity, and the 
measured drop-off of the visibility curve at short base- 
lines indicates a size of the inner rim that, when identified 
with the dust sublimation radius, yields the approximate 
value of the dust sublimation temperature. However, de- 
tailed models of the rim have indicated that this compo- 
nent by itself cannot account for the strongest observed 
NIR bumps. A second, and perhaps even more serious, 
difficulty with this interpretation is the apparent absence 
of a prominent secondary bounce in the visibility curve 
at intermediate baselines, which should be present if the 
emission region has a sharp edge. To address these two 
problems, several authors have invoked the presence of 
an extended NIR-emitting disk component of tempera- 
ture T > Tsuh within the dust subl imation radius. In 
the cas e of MWC 275, for example, iTannirkulam et ahl 
(|2008al lbl) suggested that this component most likely rep- 
resents hot gas, but iBenistv et al.l ()201Cl[ ) argued that a 
realistic gaseous disk would not be consistent with the 
spectral and interferometric constraints on this source 
and proposed that the hot inner disk component consists 
instead of refractory dust. However, the implied subli- 
mation temperature of this dust is significantly higher 
than what is commonly assumed. 

As we demonstrated in Section[31 a hydromagnetic disk 
wind, which produces a stratified density distribution of 
dusty gas beyond the dust sublimation radius, naturally 
accounts for both the spectral and the interferometric 
properties of the observed HAe systems. In particular, 
for plausible outfiow parameters, the optically thick (to 
the stellar radiation) region of such a wind subtends a 
much larger solid angle at the star than the puffed-up 
inner rim of the disk, making it possible for this model 
to account for even the strongest measured NIR excesses. 
In addition, given that the NIR-emitting region in the 
wind is generally spatially extended, the predicted visi- 
bility curves lack a pronounced secondary bounce. These 
features are in accord with the observations and contrast 
with the predictions of the spatially localized rim sce- 
nario. In the context of our simplified model, similar 
spectral fits to the NIR bump can be obtained by us- 
ing different combinations of the parameters Mout (the 
wind mass outflow rate) and Tgub (the nominal dust sub- 
limation temperature). Using the example of AB Aur, 
we showed that this degeneracy can be broken by also 
fitting the visibility curve. Although in this paper we 
focused on the application of this model to HAe stars, it 
should be relevant also to other effectively low-L, proto- 
stars, including, in particular, cTT stars. 

The hydromagnetic disk-wind model also provides a 
natural framework for interpreting the pronounced NIR 
variability exhibited by sources such as MWC 275. As 
we showed in Fi gure the NIR burst measured by 
ISitko et all (l200l can be attributed to an increase in the 
mass outflow rate (and hence in the density and optical 
depth of the wind) by ~ 75 %. In the simple self-similar 
wind model that we have adopted, this increase could be 



brought about by a commensurate increase in the mass 
inflow rate. However, in general, the disk's capacity to 
launch a centrifugally driven outflow from any given ra- 
dius is determined by several factors, which include the 
magnetic fleld-linc inclination at the disk surface (see 
Section 12. 1|) as well as the degree of ionization and the 
magnetic-to-thermal pressure ratio inside the disk (e.g., 
iSalmeron et al.|[2007[) . Each of these factors is affected 
by different physical processes that may act on different 
spatial scales, and it is therefore also possible that there 
could be pronounced changes in the local outflow rate 
that are not linearly correlated with the local behavior 
of Min- We recall in this connection our argument that 
the launching zone of the dusty outflow is effectively self- 
limiting to the narrow region between Vguh and ^ 2 r^uh 
(see Section 12.21 and Appendix |A]) . It can therefore be 
expected that any changes in the wind-launching condi- 
tions in this radially localized region of the disk would 
have a direct impact on the magnitude of the generated 
NIR excess. 

Several proto s tars, including MWC 275 
(jWisniewski et all l2008f ). have been found to ex- 
hibit photometric variability in the outer regions of 
the circumstellar disk on time scales (of a few years) 
that are much shorter than the local dynamical time. 
This behavior is most readily interpreted in terms of 
variable sh adowing by density inhomogeneities in the 
inner disk. IVinkovic fc Jurkid ()2007[ ) pointed out that 
this shadowing cannot in general be explained in the 
context of the rim model and advocated an alternative 
explanation in terms of dusty "clouds" that are ejected 
from the disk and transit across the line of sight to 
the source. The same clouds can also produce the 
observed NIR variability and, in a similar vein, may 
account for the anti-correlated near- and far-infrared 
variability that has been observed in the cTT star 
LRLL 31 ()Muzerolle et al.l [20091) . IVinkovic fc Jurkid 
(|2007| ) noted, however, that the biggest drawback to this 
picture was, in their view, "the lack of a known force 
capable of lifting a dust cloud out of the disk." This 
difficulty can in principle be overcome within the frame- 
work of the disk-wind model, where one can appeal, 
for example, to the ram pressure of the magnetically 
driven homogeneous wind component as one possible 
mechanism for uplifting clouds from the diskQ This 
process was previously considered in con nection with 
the in terpretat ion of maser disks in A GNs (jKartie et al.l 
119991 : see also iKondratko et al.l |2005() . As we noted in 
Section [l] the dusty disk-wind interpretation of the 
^ 3 /im bump was also originally proposed in the context 
of AGNs. It is interesting to observe in this connection 
that detailed studies of the infrared emission properties 
of this class of sources have provided strong evidence 
that typical disk outflows in AGNs are, in fact, clum py 
(e.g., lElitzur fc Shlosman|[200l iNenkova et al][2008l ). 

To check on the self-consistency of the scenario 
wherein a cloud is uplifted by the ram pressure of 
a disk wind, we estimate the maximum hydrogen 
column density of an uplifted cloud (subscript 'c'), 

^■^ ITambovtseva fc: GrininI I I2008I) similarly suggested that inho- 
mogeneities in the inner disk wind, which they attributed to a 
generic turbulence of the outflowing gas, could produce variable 
shadowing effects in the outer disk. 



14 



Bans and Konigl 



-^c,max = ncRc (whcrc Tic is the cloud's mean par- 
ticle density and Rc is its radius), by balancing the 
upward wind ram-prcssurc force on the cloud against 
the downward tidal gravitational force. This yields 
A^cmax «^ Mout/[(GM*r)i/2^m.p] = 2.3 x lO^^ {M,^,/5 x 
10-« Mq yr-i)(M*/2 Mo)-V2(r/0.3 AU)-!/^ cm-^ 
(where /x w 2.33 is the molecular weight and nip is 
the mass of a hydro gen nucleus; cf. Equation (24) 
in iKartje et al.l Il999[ ). A cloud with this column 
will indeed block most of the stellar radiation that it 
intercepts. We can also roughly estimate the transit 
time ttr of a cloud of this type that moves across the 
line of sight in the vicinity of the disk surface. Taking 
a typical cloud density pc and projected speed utr to 
be ~ p(r = rsuh,z = 0) and i/'o^k, respectively (see 
Equations ([3]) and Q), we obtain ^tr = ^Rc/vt^ w 1.2yr 
for our fiducial parameter values. This estimate is 
consistent with the in f erred values in a source lik e 
MWC 275 (|Sitko et al.l [200l IWisniewski eTall [2008h . 
We note, however, that alternative mechanisms for 
generating a clumpy outfiow from a magnetized disk 
— for example, a process that resembles coronal mass 
ejections in the Sun — are also conceivable. 

5. CONCLUSION 

In this paper we present a new interpretation of the 
near-infrared excess, dominated by a ~ 3 /um "bump," 
that is exhibited by the monochromatic luminosities of 
protostellar systems (notably low-mass classical T Tauri 
stars and intermediate- mass Herbig Ae and Be stars). 
The bump can be approximated by a thermal emission 
component whose temperature (^ 1500 K) is of the or- 
der of the sublimation temperature of typical interstel- 
lar grains, which has motivated the suggestion that it is 
produced in the vicinity of the dust sublimation radius 
in these systems. This interpretation is supported by 
the finding that the interferometrically measured sizes 
of low-luminosity sources scale with the bolometric lu- 
minosity as expected for Tgub- The most widely invoked 
scenario until now has been the "puffed-up inner rim" 
model, which postulates that the protostellar accretion 
disk in sources with pronounced NIR excess is optically 
thin to the stellar photons, so that the latter can prop- 
agate along the equatorial plane until they are absorbed 
by the dust in the disk at Vguh and subsequently reemit- 
ted at NIR wavelengths. In this picture, the inner edge 
of the dusty disk expands as a result of heating by the 
absorbed radiation, which increases the fraction of stellar 
photons that are intercepted and reprocessed by the disk. 
Although this model is able to relate the site of the bump 
emission to the dust sublimation radius, it has faced two 
important difficulties: It cannot account for the observed 
strength of the bump in the brightest sources, and it pre- 
dicts (on account of the inherent sharpness of the dusty 
disk's inner edge) prominent secondary peaks in the in- 
tcrferometric visibility curve, which are not observed. 

According to our proposed interpretation, the domi- 
nant contribution to the NIR bump in protostars comes 
from the reprocessing of the stellar radiation in a dusty 
hydromagnetic disk wind. In this picture, the dust is up- 
lifted from the disk by a centrifugally driven wind, and 
since dust is present only beyond rguh, this naturally ac- 
counts for the association of the bump emission region 



with that radius. Since most of the stellar photons in this 
scenario are intercepted above the disk surface, it is not 
necessary, as in the rim model, for the disk to be optically 
thin to the stellar radiation at r < rgub • Furthermore, 
since the dusty wind subtends a much larger solid angle 
at the star than the dusty disk's inner rim, this model 
can reproduce the strength of the bump even in bright 
sources such as the HAe stars AB Auriga and MWC 275. 
The lack of a pronounced secondary bounce in the visi- 
bilities measured in these sources by high-resolution in- 
terferometry is also naturally explained by the extended 
nature of the emission region in the disk-wind model. 
In fact, using a simplified wind/disk model and a spe- 
cially constructed Monte Carlo radiative transfer code, 
we showed that a single combination of physically plau- 
sible wind parameters can generate a good fit to both the 
SED and the visibility for each of the above two sources. 
In the case of AB Aur, we demonstrated how the param- 
eter degeneracy that arises when one fits the SEE) alone 
can be lifted when both the SED and the visibility are 
fitted simultaneously. 

A centrifugally driven wind along a large-scale, ordered 
magnetic field is a leading candidate for the physical 
mechanism that underlies the ubiquitous energetic out- 
fiows in protostellar systems. There are observational in- 
dications that at least some of the outflows are launched 
from the circumstellar accretion disk on scales that arc 
comparable to fgub, and that some outflows contain dust 
that likely originates in the disk. Based on the wind pa- 
rameters of our model spectral flts, we argued that the 
spatial extent of the launching region of a dusty disk out- 
flow could be self-limiting to a narrow radial range (be- 
tween Tsub and ~ 2rsub) on account of the shielding by 
the uplifted dust of stellar FUV photons that play a key 
role (through photoevaporation) in the mass loading of 
the wind. We also suggested that the disk outflow model 
provides a natural framework for interpreting the spec- 
tral variability measured in several of these sources. In 
the case of MWC 275, we showed that the reported recent 
~ 30 % increase in the 1 — 5 p,m flux could be attributed 
to a 75 % increase in the local mass outflow rate. It 
has previously been proposed that the NIR variability in 
protostars is associated with the ejection of dusty clouds, 
and we pointed out that the ram pressure force exerted 
by a homogeneous hydromagnetic wind is one plausible 
means of uplifting such clouds from the disk. 

It is interesting to note that a ~ 3 fim NIR bump is 
also observed in certain AGNs, and that the association 
between the size of the region from which this emission 
originates and the dust sublimation radius was first made 
in that context. The interpretation discussed in this pa- 
per, wherein the NIR excess is attributed to the repro- 
cessing of radiation from a central continuum source by 
a dusty disk outfiow, was also previously considered in 
connection with Seyfert galaxies and quasars. More re- 
cent modeling of the infrared emission in AGNs has es- 
tablished that the dusty winds in these sources are likely 
clumped, which brings out yet another possible similarity 
with protostellar disk outfiows. 

There is evidence that low-luminosity low- and inter- 
mediate mass protostars (comprising classical T Tauri 
stars, Herbig Ae stars, and low-L* Herbig Be stars) dif- 
fer from high-luminosity Herbig Be stars in their observa- 
tional manifestations, including the Ha line polarization 
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Fig. 8. — Spatial distribution of the stellar FUV photons that reach the disk surface, derived from a MCRT calculation for the best-fitting 
model of AB Auriga (corresponding to A/out = 4 X lO~*M0yr~^ and Tgub = 1850 K; sec Figures [3^ and|4^). The solid vertical line 
indicates the location of the adopted outer truncation radius of the disk wind, whereas the other vertical lines indicate the locations on 
the disk surface where a column density through the dusty, molecular wind to a point at a latitude of 60° on the surface of AB Aur would 
have the value marked in the figure. 



(jVink at al.ll2005[). the locatio n on the size-luminosity 
diagram ( Monnier et al.l l2005f) , and the distribution of 
visibilities (IVinkovic fc Jurkicll200l . In the context of 
the dusty disk-wind scenario, high-luminosity sources are 
distinguished by the fact that the effect of radiation pres- 
sure on grains may affect the dynamical structure, and 
consequently the radiative prop erties, of the dusty por - 
tions of their disk outflows (e.g.. iKonigl fc Kartielll99-5 ). 
We plan to investigate the extent to which this effect may 
influence the observed characteristics of high-L* proto- 



stars in future work. 
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APPENDIX 
A. DISK SHIELDING BY DUSTY WIND 

I n Sectionl2.2lwe noted that, if external heat ing of the disk surface is required for efficient mass loading of the wind 
(cf. iCasse &: Ferreiral 120001 : iPesenti et al.ll200^ . then a strong dusty outflow may shield the outer regions of the disk 
from the heating stellar X-ray and FUV radiation and could thus be self-limiting in its radial extent. For the wind and 
stellar parameters of the HAe systems considered in this paper, we suggested that the dusty outflow, which starts at 
the dust sublimation radius r^ub, would typically not extend much beyond ~ 2 rguh- Since this effect could in principle 
apply to other types of protostars and conceivably also to other astrophysical systems — such as AGNs — in which a 
disk wind, especially one that is dusty, could be self-shielding in this way, it may have important implications to our 
general understanding of astrophysical disk winds. We therefore consider this effect here in somewhat greater detail. 

The ba sic figure of merit th a t we identified is the column density to the absorption of FUV (6 — 13.6 eV) photons, 
which, as lGorti fc Hollenbachl (j2009f ). for example, pointed out, depends in detail on the opacity of the uplifted dust 
but is typically ~ lO^^cm"^. Incidentally, the absorption of X-ray and FUV photons is considered in the just-cited 
reference in connection with the photoevaporation of circumstellar disks by the stellar radiation. T his may, in fact, be a 
useful way of thinking of the effect of these photons on launching the disk outflow. As explained in lGorti &: Hollenbachl 
(j2Q09( ). the evaporative outflow that is initiated in this way would not be strong enough to overcome the downward 
pull of tidal gravity in the inner disk region that is of interest to us. We note, however, that the initiation of this 
outflow may nevertheless have an important effect on the mass loading of a hydromagnetic wind, in which magnetic 
stresses are primarily responsible for overcoming the gravitational force. 

A rough estimate of whether the dusty wind is self-shielding to FUV photons can be obtained with the help of the 
self-similar wind model that we have employed in our calculations. Considering the limiting case of a ray that grazes 
the disk surface, we can use Equations ([3]) and (|4]) to evaluate the hydrogen-nucleus column density of dusty molecular 
gas beyond Vsnb to be iVn = (2 piri//ump)(ri/rsub)"'^^^- Using n = 2.33, rsuh/i"i — 0.3, and the fiducial parameter 
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values adopted in Equation we get A'h = 5.0 x lO^^cm^^. While this is clearly an upper limit, it is useful in 
indicating that the bulk of the FUV photons that arrive at the disk from above are likely to be absorbed in the wind 
over a restricted radial range, so that the wind could be self-limiting in the sense discussed above. We now present 
the results of explicit calculations that quantify the extent of this radial range. 

In a real system, stellar photons can reach the surface of an optically thick disk on account of the finite solid angle 
that the star subtends at any point on the disk's equator, because of the flaring of the disk's surface, and through 
scattering by the dust in the wind. However, in our simplified MCRT scheme, in which the star is treated as a point 
source and the disk is taken to be flat, only dust-scattered photons reach the disk surface. We have verified, however, 
that the temperature distributions of both the dusty surface layer and the thermalized disk interior that we use to 
approximate the thermal structure of the dusty disk (sec Section [2?3l) obey the T(r) oc r~^/* relation of a fiat, optically 
thick disk that is irradiated by a finite-size star (e.g..|Adams fc S hu 1986}). This indicates that the approximation of 
a point source does not have a significant qualitative effect on the calculated thermal structure of the disk. 

Figure [8] shows the results of a MCRT calculation for model parameters appropriate to AB Auriga. The figure shows 
the spatial distribution of the FUV photons that reach the disk, which are a fraction ^ 5 — 10% of the total FUV 
photons emitted by the s tar. (For com parison, a flat, infinite, windless disk around a finite-size star intercepts 1/4 of 
the emitted photon; e.g.. lAdams &: STi u 1986,) The solid vertical line in the figure marks the location (rgmax ~ Srgub) 
of the outer wind truncation radius in our model. It is seen that this location roughly corresponds to the FWHM point 
for the distribution of stellar photons that reach the disk. The other vertical lines were obtained from a complementary 
calculation, in which the column densities through the dusty, molecular wind to different locations on the disk surface 
were calculated along rays originating at a latitude of 60° (mimicking the possible enhancement of UV emission at the 
location of a magnetospheric accretion shock on the stellar surface) on a sphere of radius 2.4 Rq (the inferred radius of 
AB Aur) centered at the origin. It is seen that the nominal outer truncation radius lies in between the disk radii that 
correspond to columns of 1.0 x 10^^ cm~^ and 1.5 x 10^^ cm~^, respectively, which is entirely consistent with the figure 
of merit for Nh adopted above. These results support the idea that the launching zone of a robust dusty outfiow in 
the disks of HAe stars could be restricted to the vicinity of r^uh by the absorption in the dusty outflow of stellar FUV 
photons that would be needed to mass- load the wind further out. 

B. MONTE CARLO RADIATIVE TRANSFER SCHEME 
B.l. Geometry 

As is standard in most MCRT calculations, the determination of events (in this case scatterings and absorptions) 
occurs in smaller subdivisions or zones of the modeled geometry. The zones are taken to be sufficiently small to ensure 
that the values of variable parameters would typically be nearly constant across them. In the present application, the 
relevant variables are the grain temperature (calculated in the code) and the gas density (which is an input variable 
obtained from the wind model). The configuration of zones generally reflects some symmetry in the problem; in the 
case of the self-similar centrifugally driven winds that we consider, the shape of the streamline (given by Equation ^) 
is the natural geometry on which the zone structure can be based. Since the density of dust in the wind decreases with 
cylindrical radius r, we expect most of the "action" to be near the inner boundary of the dusty zone at Tsub, which 
would call for a higher resolution (smaller zones) in that region. With this in mind, our zones are defined horizontally 
by streamlines launched from logarithmically sampled footpoint radii vq. We keep the vertical boundaries of the zones 
uniformly spaced. 

B.2. Initial Sampling of Source Luminosity 

Th e MCR T scheme used in this paper is based on the "packet" method described in iBiorkman fc Woo"dl (|2001[ ) and 
iLucvl ()1999[ ). The underlying principle of this method is that photons are sampled from an initial source luminosity 
Lf in the form of Np equal-energy packets. The energy of each packet is thus given by 

Ep = ^. (Bl) 

Packets are strictly monochromatic, implying that, since each packet's energy is held constant, the number of photons 
in a packet is wavelength dependent. This implies that there is a larger number of photons in a packet corresponding 
to infrared radiation than in a packet representing blue light. In the numerical scheme, the wavelength of any given 
packet is sampled by assuming that the source is a blackbody of temperature T* . Dividing the emission into K sampled 
wavelengths, equal-probability wavelengths Xk are determined by 

A- J^Bx{T,)d\ ^ ' 

where B\{T^,) is the Planck function. 

A source packet's wavelength is determined in the co de b y randomly drawing an integer k in the interval (0, K — 1) 
and then using the corresponding Afc from Equation (jB2p . The initial angle at which a packet from the source 
travels into the system (i.e., through the innermost streamline of the dusty wind region) is taken to be random and is 
determined from 

e ^ cos-\2Ci ~ I) . (B3) 
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Here, and throughout this Appendix, Q represents a random deviate in the interval (0, 1). 

B.3. Radiative Events 

As the photon packet travels through the dusty wind, there is a probability for it to interact with the dust grains 
through either scattering or absorption. Whether or not an event occurs is a function of the absorption and scattering 
opacities, and the first step in determining whether any event takes place within the zone is to calculate the effective 
maximum optical depths for these interactions: 

T^/'^ p(T,z)Kfdl (B4) 

The superscripts 'a' and 's' refer to the absorption and scattering processes, respectively. The point {tc, Zc] represents 
the current location (in cylindrical coordinates) of the packet, whereas the point {rcndj-^cnd} represent the location 
on the boundary of the next adjacent zone that a packet would hit if traveling a maximum possible path length dl 
through the current zone. The path length dl is a function of both Tc and Zc and of the packet's current direction angle 
6. [Since our scheme is formally 2D, the packet direction has no azimuthal (0) dependence.] The following condition 
determines whether there will be any event (cither scattering or absorption) within the current zone: 

C2 < 1 - 6"^'°' . (B5) 

The optical depth rtot refers to the sum of the absorption and scattering optical depths (rtot = ''"a +'''a)- ^^^^ condition 
does not hold, the packet continues traveling until it reaches {^cndj-Zcnd} and the steps outlined in Equations (|B4p 
and (|B5|) arc repeated for the new zone. If the condition (|B5|) holds and there is an event in the current zone, the 
location within the zone, {r', z'}, where the event occurs is determined by finding the effective optical depth that the 
given packet actually reaches, i.e., 

G = ^^"■'•'^ . (B6) 

Ttot 

Here Ktot is the sum of the absorption and scattering opacities. After this distance is determined, a newly drawn 
random deviate determines what type of event occurs: 

C4 < • (B7) 

_|_ 

When the condition ()B7|) holds, the event is absorption, whereas if it does not, the event is scattering. 

B.4. Scattering 

Scattering in this problem is effectively purely elastic, and in this analysis it is also treated as being isotropic. Hence 
the only variable that is changed by a scattering event is the angle at which the packet is traveling. The new angle 
after scattering is randomly drawn as in Equation (jB3p . 



B.5. Absorption 

As the code runs, the number of packets absorbed in each zone, N^^'^, is tallied and used to calculate the local grain 
temperature T^. Since each packet has the same energy, the total energy absorbed in a zone is (using Equation (|B1[) ) 

Eabs^^absj^^^^.^sLM ^ (B8) 

The total energy emitted in a zone over the time interval St is: 

E-^P = 4TTSt J dV,j J pKlBx{T,j)d\ . (B9) 

Expressing Equation (|B9|) in terms of the Planck mean opacity Kp(r) = J K\B\{T)dX/ (aBT*), where ctb is the 
Stefan-Boltzmann constant, and using the fact that / pdVij equals the mass rrnj of grains in the zone {ij}, leads to 

E°p ^ASt myKp(r„>BT^ . (BIO) 

To satisfy radiative equilibrium, the total absorbed energy in each zone must be reemitted within the same zone, and 
therefore we equate the expressions (|B8[) and (jB10[) to obtain the following formula for the dust temperature within a 
zone: 

'^ij ^ JTt • i^^^) 

Equation (|B11[) is implicit and therefore has to be solved iterativcly for Tij each time a packet is absorbed. To save 
time, the Planck mean opacities can be precalculated for a range of temperatures, and then values needed in this 
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equation can be interpolated from the precalculated files. After the packet has been absorbed, it is reemitted with 
a new wavelength. This new wavelength is determined by the temperature correction scheme of iBiorkman fc WoodI 
((2OOI . which we now summarize for the reader's convenience. 

Imagine a packet is absorbed in the zone {ij} and reemitted right away based on the local dust temperature T/^. 
When a subsequent packet is absorbed in the zone, the temperature increases to Tij = T-^ + AT. It is clear that the 
first packet was reemitted with the wrong wavelength corresponding to Tij — AT; it should have been reemitted with 
more energy because the temperature is in fact higher. The implied additional monochromatic luminosity per unit 
mass, /S.C\, is given by 

/^Cx ^Cx- C'x = nl[Bx{T,,) - Bx{T,j - AT)] , (B12) 
which, when AT /Tij is small (which can be ensured by sampling a sufficiently large number of packets), becomes 

ACx = KlAT^^^ . (B13) 

To correct for the previously emitted spectrum, the reemitted wavelength is drawn from a probability distribution that 
is based on ACx using the expression 

+ ^ Jt^^ArfT^ fc = 0,l,2....if-l (B14) 

(cf. Eauation (|B2p ). This correction is applied to every packet that is absorbed in a given zone. For computational 
efficiency. Equation (|B14p is evaluated for a range of different temperatures before the code is set to run. For each 
temperature there are K probable reemitted wavelengths, and these values are tabulated into files and then read by 
and stored in the code. This procedure yields values for the reemitted wavelength that are functions of the zone 
temperature and of a random integer fc = 0, \...K — 1 (and thus are specified by A^°(Tij)). Note that, even though the 
packet wavelength has changed, the reemitted packet still has an energy Ep, so energy is explicitly conserved. 

B.6. Spectral Energy Distribution 

Each packet that exits the system is binned into A'^Xbi,, wavelength bins and A^inc inclination bins, with indices q and 
Z, respectively: 

' = ^mcCOS0. (B15) 

Here Amax is chosen to be some maximum wavelength that the processed packets reach; for the region modeled in this 
analysis, A,nax ~ 1000 /im. The specific flux per bin of the exiting packets is 

where d is the distance to the observer and A^qi is the number of packets that occupy the {g/} bin. The total flux from 
the star, T* = L^, / [And?), can be rewritten with the help of Equation (jBl[) as 

The effective width of each bin is AA = Amax/-^Abi„i which can be ex presse d, us ing Equation (jBlsp . as A A = \/q. 
Taking the wavelength at the center of the bin and using Equation (|B16p and (|B17[) yields a formula for the normalized 
SED flux: 

^-(g + 0.5)^. (B18) 
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